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Abstract. This paper aims to validate the /3-Ginibre point process as a model 
for the distribution of base station locations in a cellular network. The /?- 
Ginibre is a repulsive point process in which repulsion is controlled by the /3 
parameter. When /3 tends to zero, the point process converges in law towards 
a Poisson point process. If /3 equals to one it becomes a Ginibre point process. 
Simulations on real data collected in Paris (France) show that base station 
locations can be fitted with a /3-Ginibre point process. Moreover we prove that 
their superposition tends to a Poisson point process as it can be seen from real 
data. Qualitative interpretations on deployment strategies are derived from 
the model fitting of the raw data. 


1. Introduction 

Statistical models of transmitters locations aim to provide tools to understand 
real network deployment. For telecommunication companies, the a priori knowl¬ 
edge of the distribution of the antenna locations helps to predict and manage the 
costs of a network deployment. Such models also provide mathematical tractable 
methods to estimate the coverage probability of a given network. These results 
would also interest telecommunication regulators and public health authorities, 
since electromagnetic exposure has become a worldwide issue. 

The first model introduced in radio networks was the regular hexagonal deter¬ 
ministic network. Although the regular lattice of cells gives an approximation of the 
cellular concept, it fails to catch the proper reality of network deployment. It proves 
also to be an optimistic bound in terms of interference estimation [T] . The random 
nature of the parameters involved in defining a proper coverage strategy makes it 
difficult to use a deterministic and regular model. Stochastic geometry ideas, espe¬ 
cially about random point processes -i.e. Poisson Point Processes (PPP), Matern 
hard-core point processes, Ginibre Point Process (GPP) and /3-Ginibre Point Pro¬ 
cess {,5-GPP)- were then widely explored in the wireless communication literature. 
Pioneer works in this field were realized by Baccelli et al. on PPP [5]. Many results 
were then derived, such as the coverage probability in respect of the SINK. Last 
developments of PPP models also include modelling of heterogeneous (fc-tier) net¬ 
works [3]. However, positions of the base stations in a PPP deployed network are 
uncorrelated with one another. Therefore clusters of points may occur. Mean inter¬ 
site distance of such configurations is thus smaller than what happens in reality. As 
a result, PPP models generate more interference than that of a real network. The 
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articles of Andrews et al. [T] and Nakata et al. [1] show that the PPP provide with 
the most pessimistic prediction of outage probability compared with other repulsive 
models. 

Spatial correlations between base stations locations exist, since they have to be 
separated from one another to maximize coverage and minimize inter-site inter¬ 
ferences. To take into account these effects, repulsive (or regular) models were 
introduced in the literature. A simple approach is to transform a PPP into a repul¬ 
sive point process by thinning. Such processes are called Matern hard-core point 
processes. Interferences for such deployed networks were investigated in but 
hard-core models proved to be difficult to manipulate since the outage probability 
can not be analytically deduced. Soft-core processes then rose community’s in¬ 
terest. Among them, GPP and /3-GPP (two determinantal point processes) were 
investigated in the wireless communication field. They were at first introduced 
by Shirai et al. [6] in quantum physics to model fermion interactions. Works of 
Miyoshi et al. [7] and Deng et al. [5] have derived coverage probability in respect of 
the SINK for both GPP and /3-GPP models. 

In this paper, we show that base station distribution for an operator and for 
a technology can be fitted with a /3-GPP distribution in the Paris area. The dis¬ 
tribution of all base stations of all operators can be fitted with a PPP. Our main 
contribution lies in the theoretical justification of this phenomenon. We prove that 
the superposition of different /3-GPPs converges in distribution to a PPP process. 
Finally we draw conclusions on the coverage-capacity trade-off made by different 
operators. Qualitative results are derived from the inferred values of /3 and the 
intensity A. A can give information on the dimensioning strategy adopted by the 
operator, while /3 give insights on the coverage. 

Other existing papers on antenna deployment models mainly consider the com¬ 
putation of the SINK and coverage probability for a wide set of point processes. We 
are instead interested in validating the /3-GPP model and the PPP superposition 
model with real data on a dense urban area. Such a case study is made possible 
because French frequency regulator (ANFR) provides location in an open access 
database [9]. 

In Section II, we define mathematically the /3-GPP and introduce the conver¬ 
gence in distribution theorem for a superposition of /3-GPP. In Section III, we give 
the method used to fit the /3-GPP model with the actual data. A qualitative in¬ 
terpretation of the deployment strategies is then realized from inferred /3 and A 


2. Theoretical Model 

In this section, we recall the definition of the ,5-GPP. We also introduce the 
theorem /3-GPP convergence theorem. 


2.1. Definitions. Let C denote the complex plane. Let $ be a realisation of a 
certain point process. Let xi,... ,Xk, be fc-tuples of distinct pairwise elements of 
C. Let i? C C be a Borel set and / : i?" —>■ be any Borel function. 


Definition 1 (Gorrelation function). The fc-th joint density function of a point 
process is defined by: 


E ^ f{xi,...,Xk) 

(£Ci,...a:„) 

Xi^Xj 


/ f {xi,... ,Xk)p’'’'Hxi,.. .,Xk)dxi.. .dxk- 
Jw' 
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Let /3 be a real number in ]0,1], let A, a strictly positive real number be the 
intensity of a point process and c = Att. The /3-GPP is a determinantal point 
process that can be defined by its correlation functions. 

Definition 2 (/3-Ginibre). 

{xi, ...,Xk)= det (Lfc./9 {xi,Xj) ,1 <hj < k), 

where Kc^p is a kernel such as V(a;, y) G C^: 

Kc,p {x,y) = 

TT 

with the respect of the Lebesgue measure. 

A /3-GPP can also be obtained by thinning and rescaling a GPP. Each point 
of the GPP is kept independently with probability /3, then a rescaling of ratio 
is applied in order to maintain the original intensity. If /3 = 1, the 1-GPP is 
equal to the GPP. For /3 tending to zero, the thinning increases the randomness 
of the process, and then the kernel tends to a diagonal matrix. Gorrelation 
between points disappears. The Laplace transform of a /3-GPP is given for all / 
by L{f) = det(/ — Vl — e~fKc^pVl — e“-f), where the determinant is the one of 
Garleman-Fredholm. It is easy to see that the /3-GPP Laplace transform converges 
then to the one of a PPP when /3 goes to zero. 

2.2. Properties of superposition of /3-GPP. One of the main novelties in this 
paper is the study of the superposition of multiple /3-GPP. We give the key conver¬ 
gence theorem for the /3-GPP. 

Let $ 1 , $ 2 , be point processes in a metric space R. Let y be the space of all 
locally finite subsets (configurations) in R. We now introduce the /3-GPP superpo¬ 
sition convergence theorem. 

For all n G N*, let be the superposition of ^|j family of /3„,i-GPP with 

intensity Xn,i = {mT) ^ Ci and Pnp G]0, 1]. 

Theorem 1 (/3-GPP superposition convergence theorem). Let us suppose that: 

(i) the sequence (ci)igN* C is bounded, 

(ii) lim n~^ Ci is finite and equal to c. 

n—>-+cx) 

Then ($rt)„gisj. converges in distribution to a PPP $ with intensity C7r“^. 

Proof. This theorem is proven in Appendix |A] □ 

Hypotheses (i) and (ii) of ThU] are quite restrictive because the intensities of 
each /3-GPP are dependent of n. However, in practice, we mainly work with finite 
families of /3-GPP. Therefore, we can choose the value of the (ci)igi...„ such that 
they match the real values of the intensity of each /3-GPP. 

3. Simulations 

In this section we introduce the fitting method that is used to obtain the pa¬ 
rameter /3. We also present the results from the fitting of each deployment and 
operator in Paris, France. 

3.1. Summary statistic. In order to fit the real deployment to the /3-GPP model, 
we introduce the J function that characterize any point process. This function is 
a summary statistic based on inter-point distances. General information about 
summary statistics can be found in m- Let u be any location in the plane C, and 
$ = {xij-gpj be a realization of a /3-GPP. 
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Definition 3 (J-function of the /3-GPP). 


Vr e k; 


J{r) 


1 - G{r) 

■ 1 - F{r) ’ 


where G and F are respectively the contact distribution function and the nearest- 
neighbor distance distribution function. 


The J-function characterizes the repulsiveness or attractiveness of a point pro¬ 
cess. If the J-function is bigger than one, the point process is repulsive, otherwise, 
it is attractive. By definition of a PPP, J(r) = 1 In the case of the /3-GPP, the J- 
function is always bigger than one on its definition domain. This confirms that the 
/3-GPP is a repulsive point process. It characterizes the ,S-GPP and its expression 
is proved in [5]. 



Figure 1. Example of data sample for one GSM operator. The 
J-function is fitted on the points within the rectangular window. 

3.2. Fitting method. Thanks to R language and the spatstat package the 
estimate of the J-function is derived from the raw data. Since we consider only 
a finite set of antennas, edge-effect might appear on the J-function estimate. We 
then have to keep a subset of the data to perform the estimation. Figured] gives the 
window we considered for extracting data in Paris, France. It covers about 60% of 
the city and its shape is chosen to match the geographical borders. The values of 
the J-function estimate are computed for r < 600 m. Above 600 m, the estimation 
is not relevant due to the edge-effect. J is then directly fitted on the estimate and 
the parameter /3 is deduced. An example of fitting is given in Figured) It is clear 
that the point process formed by the base stations locations is repulsive and fits well 
the theoretical model. Therefore, it outfits the PPP model, because the J-function 
a PPP is equal to one for all r. In the next paragraph we present the results we 
obtained on raw data. 

3.3. Fitting results and interpretation. Locations of the base stations are pub¬ 
licly available for the whole French territory and can be found online [5]. There 
are four operators in France and most of them provide 2G to 4G coverage. For 
each operator and each technology, numerical values of P and A from the fitting are 
given in Table [T| 

Values of /3 and A give some insights about the deployment strategy of each cel¬ 
lular network operators, especially about the coverage-capacity trade-off. Orange’s 
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Figure 2. Example of J-function fitting for Orange, SFR and 
Bouygues on the 3G 2100 MHz band. As a comparison, J(r) = 1 
for all r in the PPP case. 


Table 1. Numerical values of /3 and A per technology and operator 



Orange SFR Bouygues Free 

GSM 900 ^ 

A 

0.81 0.76 0.65 NA 

2.39 2.65 2.63 NA 

GSM 1800 ^ 

A 

0.84 0.85 0.71 NA 

3.00 2.39 3.59 NA 

UMTS 900 ^ 

A 

NA 0.97 0.53 0.89 

NA 1.92 2.44 1.05 

UMTS 2100 ^ 

A 

1.04 0.65 0.82 0.89 

3.27 3.48 4.04 1.05 

LTE 800 ^ 

A 

1.02 0.93 0.67 NA 

0.67 1.65 1.87 NA 

P 

LTE 1800 

A 

NA NA 0.75 NA 

NA NA 3.46 NA 

13 

LTE 2600 

A 

0.93 0.67 0.63 0.89 

2.80 2.76 2.46 1.05 


high values of /3 and A suggest that this operator deployed (as the historic, pre¬ 
viously state-owned operator) a network that fulfilled an optimal coverage and an 
optimal traffic capacity (densely deployed network). However, SFR and Bouygues 
first deployed a network with a minimum of antennas (in order to abide by the cov¬ 
erage requirement of the regulator) and then gradually increased traffic capacity 
on hot-spots (by increasing locally the number of antennas). This involves adding 
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more antennas on sites that are already covered, thus creating clusters and decreas¬ 
ing the value of j3 and increasing the value of A. The French telecommunication 
regulator (ARCEP) published yearly reports [12] that suggest such evolution. 

We deduce that French operators used two different deployment strategies. The 
first strategy consists in fulfilling both coverage and optimal traffic capacity at 
once. While the second strategy is to deploy a network that abides to the coverage 
requirements in a first stage, then in a second stage to increase the number of 
antennas on hot-spots in order to improve the traffic capacity. 

Table 2. Numerical values of /3 and A per operator and for the 
superposition of all the sites 



Orange 

SFR 

Bouygues 

Free 

Superposition 

p 

0.94 

0.70 

0.81 

0.89 

0.17 

A 

3.48 

3.70 

4.23 

1.05 

10.28 

Number of sites 

185 

197 

225 

56 

547 


When deploying their 3G or 4G networks, operators reused and shared some 
existing 2G sites. Therefore, we consider that classifying the base station sites per 
operator is more relevant than classifying them by technologies. Table [5] summaries 
these results. As expected, previous conclusions still hold as values of /3 are stable 
between the two tables. We also notice that Free, as a newcomer (2012), has a small 
amount of traffic to deal with, and therefore has deployed less antennas than its 
competitors. Data analysis also shows that the superposition of all sites is tending 
to a PPP as P is equal to 0.17. Therefore the PPP model still holds as an indicator 
of electromagnetic exposure of cellular networks. 

4. Conclusion 

In this paper, we successfully show that /3-GPP is a realistic model for base 
station distribution. The /3 parameter is inferred by using statistical tools on real 
data. Qualitative results on network deployment are then derived. We also prove 
theoretically that the superposition of multiple /3-GPPs converges in distribution 
to a PPP justifying observations made on real deployments. This will have greater 
implications in modelling multi-tiers networks. We show that the values of A and /3 
are characteristics of the coverage-capacity trade-off. Future works will investigate 
the impact A and /3 on the design of optimal deployment strategies. 

Appendix A. Proof of Theorem 1 

Let A be a compact subset in C. For a realization of a point process 4>, the 
random variable $(A) is the number of points in the compact A. 

Theorem 2. [Convergence in distribution theorem] For any A compact subset in 
C, if the three following properties hold: 

(i) lim P(4>„(A) = 0) = P($(A) = 0) 

n—>-+oo 

(ii) limsup P(<1>„(A) < 1) > P($(A) < 1) 

n—¥-\-oc 

(iii) lim limsup P(<1>„(A) > t) = 0 

i—H-oo 
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Thin is achieved if all conditions of Thf^] are satisfied. 

Condition (in). Thanks to Markov inequality, and since is bounded, 

(iii) holds. 


Conditions (i) and (ii). For a Poisson point process, we know that: 

P($(y4)=0) = e“l^l^’"“\ 

P($(7l)<l) = e-l^l^’^”'(l + |A|c7r-i) . 

We have yet to calculate the left-hand side of both inequalities (i) and (ii). Let Kn^i 
be the kernel of a /3„^i-GPP. Proposition 3 of Goldman’s paper [T3] states that: 

= = 1 + f det[K^,,]{vu:;Vpyndv), 

P>1 ^ AP 

P($„,,(kl) = l)=P($„,,(kl)=0) J RnA^yidz), 

A 

where £ designates the Lebesgue measure and 

Rn,iiz) = Kn,iiz, z)+y^^ K^diz, z). 

i>2 

By hypothesis of Thin (cO^gpj. is bounded. We also know that ||ifn,i||oo = Ci{mT)~^. 
We can then prove recursively for all p > 1, there exists a M > 0 such that for all 
i € N*, 

/ mY 

0 < det[iL„y](wi, ...,Up) < \\Kn,i\\^ < ( — ) . 

\mT / 

Therefore there exists two bounded sequences (e„)„gN* and independent 

of i and , such that: 

P($„,,(A) = 0) = l-^+n-"e„., 

TT-TT 

P(d>„,,(kl) = l) = ^+n-2e;. 

riTT 

Hence, 

P($„(A) = 0) = ^, 

P(<i>„(kl) = 1) > ^ ^ + a(l). 

riTT 

Therefore, 

lim P($„(y4) = 0) = 

n—¥oo 

limsupP($„(A) = l) > 

n—>-oo 'TT 

consequently (i) and (ii) hold. □ 
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